Ground-state energy and stability limit of small 3 He drops 
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Small and stable drops of 3 He atoms can only exist above a minimum number of particles, due 
to the combination of the 3 He atom Fermi statistics and its light mass. An accurate estimation of 
this minimum number using microscopic theory has been difficult due to the inhomogeneous and 
fermionic nature of these systems. We present a diffusion Monte Carlo calculation of 3 He drops 
with sizes near the minimum in order to determine the stability threshold. The results show that 
the minimum self-bound drop is formed by N = 30 atoms with preferred orbitals for open shells 
corresponding to maximum value of the spin. 

PACS numbers: 67.55.-s,36.40.-c,02.70.Ss 
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Liquid 3 He drops offer a unique combination in nature 
of Fermi statistics, neutral charge and self-bound char- 
acter [ij, |2j ■ The interest in their physical knowledge ex- 
plains the continued effort from both the experimental |3J 
and theoretical sides 0, la, la S @ towards a better de- 
scription of their particular properties. Experimentally, 
3 He drops are currently generated in the laboratory by 
means of a free jet expansion from a stagnation source 
chamber through a thin walled nozzle [1| . The estimated 
temperature of the Fermi drops is 0.15 K and therefore 
they are in their normal state. The non-superfluid char- 
acter of 3 He drops has been detected in a series of ex- 
periments [9j where single molecules were embedded in 
the drops. If the molecule is surrounded by 4 He, the ro- 
tational spectrum presents a sharp structure which has 
been attributed to the superfluid nature of 4 He, whereas 
for a 3 He drop a broad peak is observed. 

The smaller mass of the 3 He atom, and more funda- 
mentally its Fermi statistics, introduces sizeable differ- 
ences in the stability of the drop with respect to 4 He. 
Two 4 He atoms form a bound state but a minimum num- 
ber of 3 He atoms is necessary to form a self-bound sys- 
tem |5j. This difference is observed in the experimen- 
tal setup since very small 4 He drops are detected in the 
jet, starting from the dimer, whereas it has been proven 
to be difficult to observe 3 He drops with less than 1000 
atoms J3j. Reducing this number and generating drops 
with a number of atoms closer to the threshold limit is 
one of the most important challenges for the next fu- 
ture. Small 3 He drops are expected to present "magic 
numbers" [4| whose experimental determination would 
be a clear signature of their Fermi statistics and a new 
benchmark for quantum many-body theories for inhomo- 
geneous systems. 

As a result of the complexity arising from the combi- 
nation of Fermi statistics and inhomogeneity, the num- 
ber of microscopic works on 3 He drops is significantly 
smaller than those devoted to 4 He drops. The first 
systematic study was carried out by Pandharipande et 
al. |5j in the eighties using the variational Monte Carlo 



(VMC) method. This calculation used a trial wave func- 
tion incorporating backflow correlations to correct the 
nodal surface of the noninteracting system and predicted 
that a drop with 40 atoms is self bound. More re- 
cently, Guardiola and Navarro [3 carried out a detailed 
VMC calculation of small 3 He drops including in the 
trial wave function configuration-interaction-like corre- 
lations. These new correlations improved the energy in 
a significant way and the smallest bound drop was esti- 
mated to be the one formed by 35 atoms. Recently, the 
same authors |8j have obtained a lower value (32) for this 
upper-bound threshold in a diffusion Monte Carlo calcu- 
lation (DMC) restricted to the atom's number range 31- 
34. Therefore, the improvement of both the trial wave 
function and the theoretical approach has progressively 
reduced the minimum number of atoms required for a 
self-bound drop. It is worth mentioning that using a non- 
local density-functional approach |£j] this critical number 
was estimated to be slightly smaller (29) than these mi- 
croscopic calculations. 

In this work, we present a DMC calculation of small 
3 He drops around the threshold limit for self-binding. In 
the simulation we use the fixed-node (FN) approxima- 
tion [lOJ , which provides an upper bound to the exact 
eigenvalue, and the release-node (RN) method [llj to 
estimate the quality of the FN upper-bound. The ap- 
proach is the same we followed in the past in the DMC 
calculation of the equation of state of bulk [12J and two- 
dimensional |l3( 3 He. The results show that the min- 
imum number for a self-bound drop is 30 and that in 
open-shell configurations the optimal energy corresponds 
to maximum spin. The latter conclusion was also pointed 
out in previous density functional ||| and VMC work |j| • 

The sign problem in a DMC simulation is usually 
dealt within the FN approximation |lfj. Along the cal- 
culation, the wave function /(R, r) = V't(R)^'(R, t) 
(R = ri,...,rjv) evolves according to the imaginary- 
time (r) Schrodinger equation, with ifa acting as impor- 
tance sampling function and nodal constraint. For long 
enough time r — > 00, ^(R) is the lowest energy state 



compatible with the nodal surface imposed by ^t(R)- 
The trial wave function we have used has a Jastrow-Slater 
form, 
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with a two-body correlation factor, 
f(r) = exp 
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In Eq. (JTJ, $| ($|) is a Slater determinant for the spin- 
up (spin-down) particles filled up with single-particle or- 
bitals corresponding to the polynomial part of the har- 
monic oscillator basis. From a practical point of view, 
the use of this basis is clearly advantageous since the re- 
sulting Slater determinant is of the Vandermonde type 
and therefore translationally invariant. With this model, 
and considering the Jastrow part as a function of relative 
distances only, spurious energy contributions due to the 
movement of the center of mass are not present. The 
orbitals are chosen in its Cartesian coordinate represen- 
tation, and for incomplete shells we have followed the 
prescription used by Guardiola and Navarro [2J which 
warrants invariance under 90° rotations with respect to 
the Cartesian axis. 

The quality of the upper bound in the FN approach 
depends on the accuracy of the nodal surface defined 
by V't(R-)- The nodal surface of ^x(R) © corre- 
sponds to the one of a non-interacting system. There- 
fore, this model is expected to be too simple for describ- 
ing a correlated liquid as 3 He. Using the imaginary-time 
Schrodinger equation, one can prove that the first correc- 
tion (corresponding to a short imaginary-time interval) 
to the non-interacting nodal surface corresponds to a dis- 
placement of the coordinates due to correlations with all 
the other particles [15J. These corrections are known as 
backflow correlations and are constructed by replacing 
the coordinates r.; of particles in the Slater determinants 
by 
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Similarl y t o previous studies of the homogeneous liq- 
uid [l2, llj|, we have used for the function r)(r) a Gaus- 
sian, r/(r) — exp[— ((r — tbVwb) 2 ]; Ab, ?~b, and wb are 
variational parameters. 

The FN method provides a rigorous upper-bound on 
the ground-state energy but does not provide information 
on the quality of the upper bound, i.e., the difference be- 
tween the energy obtained and the exact eigenvalue. In 
order to estimate the bias due to the particular model 
nodal surface we have used the released node (RN) tech- 
nique [ill • I n the RN method, the walkers are allowed to 
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TABLE I: Total (E) and kinetic (K) energy per particle of 
small 3 He drops as a function of the number of atoms. S z 
is the z component of the total spin of the drop. Figures in 
parenthesis are the statistical errors. 



cross the nodal surface determined by i/'t (R-) for a finite 
lifetime t r and a sign + or — is assigned to each one. 
To this end, the importance sampling wave function is 
chosen positive 



V(R) = (^t(R) 
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with a a constant, and the fermionic energy is obtained by 
projecting on the antisymmetric component i/jt(R). The 
method would arrive to the exact ground-state energy for 
t T — > oo, but this limit is not accesible in liquid 3 He due 
to the rapid emergence of bosonic noise. Nevertheless, 
the initial slope can be well determined and its value can 
be used for comparing different nodal surfaces and for an 
estimate of the magnitude of the bias introduced by the 
FN approximation 12]. 

All the DMC simulations have been carried out using 
the HFD-B(HE) Aziz potential |14j. which has proved 
high accuracy in the microscopic description of the bulk 
phases of liquid 4 He and 3 He 15]. The parameters enter- 
ing into -0t(R.), Eqs. J2J) and ©, are adjusted variation- 
ally. The dependence on the number of atoms in the drop 
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FIG. 1: Total energy of 3 He drops as a function of the number 
of atoms -/V. The error bars are smaller than the symbol size. 
The line on top of the DMC data is a guide to the eye. 



is only significant for the parameter a (0): for N = 30, 
a = 4.0 cr" 1 and for N — 40, a — 4.4 er -1 , increasing 
linearly with N (a = 2.556 A). The optimal values of 
the rest of parameters are: /3 = 1, v = 5, b = 1.14 <r, 
A B = 0.34, r B = 0.75 er, and uj b = 0.54 a. 

Table I contains results for the total (E/N) and kinetic 
energy (K/N) per particle as a function of the number of 
atoms N in the drop. With the exception of the first row 
for N — 32, 35, and 38 the calculations have been made 
using orbitals with invariance under 90° rotations with 
respect to the coordinate axis. According to our results, 
the threshold limit for a self-bound drop is N = 30 and 
all cases studied with equal N but diferent spin S z show 
a preferred state corresponding to the maximum value of 
the spin. The latter result can be also interpreted tak- 
ing into account the magic numbers which close a shell 
(of the spin- up or spin-down atoms) , which in the range 
studied correspond to values 10, 20, and 35. The re- 
sults contained in Table I show that the optimal energies 
follow the rule of having at least one closed shell, with 
preference for the smallest one: 10 for N — 29 and 20 for 
AT = 31 -39. 

The lowest total energies of the 3 He drops are shown 
in Fig. 1 as a function of their number of atoms N. The 
line on top of the DMC data correspond to polynomial 
fits and are only intended as to guide the eye. As one 
can see in the figure, the behavior of the energy with 
N is not monotonous in the regime studied showing a 
kink for N — 40, a doubly magic number JV-j- = JVj. = 
20. This kink is a remnant of the shell model chosen 
to describe the antisymmetry of the system in absence 
of correlations. Dynamical correlations induced by the 
interatomic potential smooth significantly this effect but, 
for these small drops, it is still clearly observable. On the 
other hand, Fig. 1 shows that in the regime N = 30 — 40 
the DMC data display a regular behavior, which is well 
reproduced by a second-degree polynomial. This feature 
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FIG. 2: Evolution of the energy with released time t t for a 
N = 40 3 He drop. Circles and squares correspond to Jastrow- 
Slater wave functions with and without backflow correlations, 
respectively. In both cases, we subtract to the RN energies 
the energy at initial time Eto = E(t T — 0) for an easier com- 
parison. The statistical error bars are essentially constant 
in this t r range; they are not shown for major clarity of the 
released signal. 



yields us confidence on the calculation itself and on the 
method followed for the selection of orbitals in the Slater 
determinants. 

Although the search for an exact and stable quan- 
tum Monte Carlo algorithm for solving the TV-fermion 
Schrodinger equation continues, the intrinsic difficulty of 
the problem raises the question about the maximum in- 
formation one can obtain at present from the available 
Monte Carlo methods. The only really stable method, 
which can manage a significant number of fermions, is 
FN. With FN one is able of computing rigorous upper 
bounds to the ground-state energy, with the only con- 
straint of the model nodal surface contained in ifa- The 
introduction of backflow correlations in the model has 
proven to be of crucial importance in order to decrease 
the bias introduced by ifa. Also for 3 He drops, this intro- 
duction allows for a much better description: for N = 40, 
backflow correlations make the energy per particle to de- 
crease 0.03 K, roughly a relative improvement of 25%. 

As in a previous work on bulk liquid 3 He, we have used 
the RN technique to get some insight on the quality of 
the upper bound provided by the FN method. Although 
the RN method is unstable in the sense that the subja- 
cent bosonic component asymptotically overwhelms the 
Fermi signal, the initial slope of the energy with the re- 
leased time i r can be determined with enough precision 
to make a com par ison between different model nodal sur- 
faces possible [l2| . This kind of analysis is shown in Fig. 
2, where the RN energies calculated with a V't containing 
or not backflow correlations are plotted as a function of 
t T - The Figure shows a reduction of the slope by a fac- 
tor of two when backflow correlations are present in ip^ 




FIG. 3: Threshold limit for self-binding in He drops. Solid 
circles correspond to the present results. Open circles stand 
for the FN-DMC calculation from Ref. |jj], and squares for 
the VMC results from Ref. Q. 



pointing to a significant improvement of the upper bound 
(an exact wave function would show zero slope). How- 
ever, this improvement is worse than the one we observed 
in the past in a similar calculation of bulk 3 He |l2j. In 
fact, also the variance observed in the present simulation 
of the drops is larger than the one estimated in bulk. 
Both comparisons point to a somehow incomplete treat- 
ment of the inhomogeneity inherent to drops. Possible 
improvements on this line could be obtained by changing 
the functional dependence of correlations, from (|ri — rj|) 
to (i"i,rj) |l6j |. Notwithstanding, this modification would 
introduce additional complexity in the calculation and 
require from new correlation functions which at present 
are not very well known. 

The present DMC results for 3 He drops, near the 
threshold limit for self-binding, are compared with re- 
cent MC data in Fig. 3. In the Figure, the VMC en- 
ergies correspond to the best variational calculation up 
to date 13 . The high quality of the upper bounds there 
achieved comes from the introduction of configuration- 
interaction correlations in the Jastrow part. The result- 
ing limit for self-binding was N — 35, quite close to our 
result (N = 30), and the difference with respect to the 
present results increases slightly with N. The orbitals 
we have used, which are essentially the same used in this 
VMC estimation, do not have in general good angular 
momentum quantum numbers. In principle, that can be 
considered a defficiency of the model and the results ob- 
tained could be worse than the ones obtained with better 
wave functions. However, Guardiola |l7| proved using 
VMC that the energies of 3 He drops are independent of 
the orbital angular moment and that they depend only 
on the spin. 

Recently, Guardiola and Navarro y| have reported FN- 
DMC results for 3 He drops in the range N = 31 — 34 (also 
shown in Fig. 3). They use the same interatomic po- 



tential and the same Jastrow factor but a different shell 
structure in the Slater determinant, and also a different 
form for the backflow function. Their results show a min- 
imal number N — 32 for a bound drop and energies that 
are above our results in the number range considered. 
We have verified that a significant part of the difference 
between that calculation and the present one lies in the 
different functional form used for backflow correlations. 
In Ref. ; 8], the form rj(r) — A/r 3 instead of a Gaussian 
is used; for a N — 34 drop it supposes an energy loss of 
-0.26 K. 

In conclusion, we have carried out an accurate calcu- 
lation of small 3 He drops using FN- and RN-DMC that 
has allowed for a reduction in the threshold limit for self- 
binding; our results show that this minimum number is 
N = 30. Similarly to previous density-functional |6J and 
microscopic calculations p| , our results confirm that the 
ground-state energy is achieved for maximum spin or, in 
other terms, for shell configurations where the number 
of atoms of one of the two species (up or down) closes 
a shell. In the N range studied, the energy shows kinks 
for N — 30 and N — 40 which correspond to magic num- 
bers of the underlying shell model; its signal is however 
quite depressed by the relevance of 3 He- 3 He dynamical 
correlations. 
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